Stereoscopic electron spectroscopy of solar hard X-ray flares with 

a single spacecraft 
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ABSTRACT 

Hard X-ray (HXR) spectroscopy is the most direct method of diagnosing energetic electrons 
in solar flares. Here we present a technique which allows us to use a single HXR spectrum to 
determine an effectively stereoscopic electron energy distribution. Considering the Sun's sur- 
face to act as a 'Compton mirror' allows us to look at emitting electrons also from behind the 
source, providing vital information on downward-propagating particles. Using this technique we 
determine simultaneously the electron spectra of downward and upward directed electrons for 
two solar flares observed by the Ramaty High Energy Solar Spectroscopic Imager (RHESSI). 
The results reveal surprisingly near-isotropic electron distributions, which contrast strongly with 
the expectations from the standard model which invokes strong downward beaming, including 
collisional thick-target model. 

Subject headings: Sun: flares — Sun: particle emission — Sun: X-rays, gamma rays 



1. Introduction 

Energetic electrons and ions have long been 
considered (Ellison and Hoyle, 1947) as possibly 
playing a key role in energy release by magnetic 
reconnection, in phenomena such as solar flares 
(Aschwanden, 2002), which are archetypes of mag- 
netic explosions and particle acceleration in the 
cosmos, as well as being central to space weather 
and its terrestrial influence. 

The spectrum of energetic electrons and the 
anisotropy of source electrons are normally in- 
ferred separately. Hard X-ray (HXR) directivity 
itself has, until now, been measured using two 
different techniques. The first is direct simulta- 
neous measurement of flux spectra from at least 
two spacecraft at well separated locations (Li et 
al, 1994; Kane et al 1988). The second method is 
based on the statistical study of the distribution 
of HXR fluxes and spectral indices over heliocen- 
tric angle (Datlowe et al. 1977; Bogovalov et al. 
1985; Vestrand 1987). However, neither of these 
methods provides direct information about down- 
ward emitted photons which are crucial in model 



testing. 

High resolution HXR spectrometry with RHESSI 
enables reconstruction of volume-averaged source 
electron flux spectra (electrons cm~^ s~^ keV~^ ) 
from bremsstrahlung HXR emission spectra (pho- 



tons electrons cm 



keV ) observed at the 



Earth (Plana et al. 2003; Kontar et al. 2004). 
However, the observed spectrum is contaminated 
by an albedo component, due to Compton back- 
scattering ('reflectivity') in the dense photosphere, 
of those primary bremsstrahlung photons which 
were emitted downward (Tomblin, 1972; Bai and 
Ramaty, 1978). Until now this albedo contribution 
has been regarded as a nuisance to be corrected 
for, either following Bai and Ramaty (Johns and 
Lin 1992, Alexander and Brown 2002;), or self- 
consistently (Kontar et al, 2006), before spectral 
inference of electron spectra. 

In this Letter we show that, on the contrary, the 
albedo spectral 'contaminant' in fact offers very 
valuable insight into the anisotropy of the flare 
fast electron distribution. It does so by provid- 
ing a view of the HXR flare from behind, like a 
dentist's mirror, except that the solar albedo mir- 
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ror is spectrally distorting so its contribution to 
the overall spectrum can be distinguished. This 
reflectivity has a broad spectral peak in the SO- 
SO keV range (Bai and Ramaty, 1978; Kontar et 
al, 2006). It decreases fast at low energies due to 
photoelectric absorption while high energy pho- 
tons are more likely to be lost to an observer be- 
cause they penetrate too deep into the solar at- 
mosphere to be scattered back to an external ob- 
server. On the other hand, an important prop- 
erty of bremsstrahlung emission is that it is al- 
ways monotonically decreasing for any electron 
spectrum, typically having a roughly power-law 
shape of spectral index which is nowhere less than 
1 (Koch and Motz, 1959). Thus, the primary emis- 
sion and the albedo 'bump' spectral signatures are 
very distinct. The observed spectrum in the ob- 
server's direction should contain an albedo bump 
feature, the strength of which is an indicator of the 
degree of downward beaming of the electron distri- 
bution. By use of this solar 'mirror' we can achieve 
a degree of 2-directional beam electron beam spec- 
trometry from single spacecraft photon spectrom- 
etry. 

We show furthermore that this insight con- 
strains the directivity of flare electrons, so strongly 
that the conventional solar flare models with 
downward beaming are excluded. Note that the 
inferred electron spectra and distributions cannot 
be explained by the coUisional scattering (Trub- 
nikov 1965) as assumed in coUisional thick-target 
(CTT) model (e.g. Brown 1971; Brown 1972). 
In the CTT scenario, even allowing for coUisional 
scattering, the downward-propagating electrons 
emit HXR bremsstrahlung which is also substan- 
tially downward-coUimated. 

2. Photospheric albedo as a natural con- 
straint on downward beaming of elec- 
trons 

For an inhomogcncous bremsstrahlung source 
of volume V, plasma density nlr), and electron 
flux spectrum F{E,S},',r) per unit solid angle in 
direction Jl' can be averaged over the source vol- 
ume. The photon flux spectrum at Earth distance 
R in direction O is 

— p p oo 

^W = n?j/ / F{E,n')Q{n,n',e,E)dEdn' (1) 



where Q{n,ft',e,E) is the bremsstrahlung cross 
section differential in e and H; n = J n(r) dV 
is a volume averaged density, the density- weighted 
mean electron flux is F{E, f2') = {nV)-'^ J n{r)F{E, n' , r) dV, 
and e, and E are photon and electron energies cor- 
respondingly. In fact Q depends only on the angle 
n'n between the incoming electron O' and the 
emitted photon O directions. 

Equation ([1]) can be approximated using a two 
directional representation based on the fact that 
here we are concerned with emission upward (u) 
toward the observer and downward (d) toward the 
scattering photosphere. The flux toward the ob- 
server /o(e, d) from a source with heliocentric angle 
9 can be written 

Io{e) = ^ nV j [0^(6, E)Fu{E) + 

' Q''{e,E)Fa{E)]dE, (2) 

where Fu^d = {nV)'^ J Fu,diE,r) n{r)dV and, us- 
ing axial symmetry, we have introduced 

^^^'^'^°^ " cos(^o - A0) - cos(0o + AO) 

/ Q{e,E,e')sm{0')d0' (3) 

the cross-section averaged over [6*0 — A6', 6*0 + A9] 
and centered at angle 9o- Hence, Q^{e,E) = 
Qie,E,9o = 0) and Q^{e,E) = Q{e,E,9o = 
180° — 9), where 9 is the heliocentric angle of the 
source. Electron spectrum F{E, 9) is defined in a 
similar two directional approximation: Fu{E) and 
Fd(E) are the density weighted volumetric mean 
flux spectra of electrons directed towards the ob- 
server 'upward' and downward respectively, also 
averaged over A0. 

The probability of photon emission along the 
direction of motion is around ten times higher 
than in the opposite direction and four times that 
at right angle, the cross-section having a typical 
angular scale A.9 ^ 30 — 50° (Koch and Motz, 
1959) for energies 50-300 keV. Moreover, the elec- 
tron distribution should have some additional an- 
gular spread (e.g. because of collisions (Trub- 
nikov, 1965; Brown 1972, Leach and Petrosian 
1983, MacKinnon and Craig, 1992) and magnetic 
fleld convergence). Therefore, for the spatially 
averaged distribution it is natural to deal with 
angular distribution averaged over A0. In our 
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calculations we take — 45° as a character- 
istic scale of angular dependency. (In the limit 
= 180° the cross-section becomes angle aver- 
aged and equation ^ is reduced to the simplified 
equation for isotropic electron distributions used 
in, e.g., Brown et al., (2003)). 

X-ray photons propagating downwards undergo 
Compton backscattering and absorption in the 
dense photosphere, the spatially integrated re- 
flected flux being expressible via a Green's func- 
tion approach (Magdziarz and Zdziarski, 1995; 
Kontar et al, 2006). 

/oo 
G(M,e,e')^d(e')'^e' (4) 

where /(i(e) is the the downward directed flux and 
G{^,e,e') is the angular (^ = cos(6')) dependent 
Green's functions (Kontar et al, 2006). The re- 
flected flux can be written 

/.(e,/i)=4^ i G(M,e,0&' 
[Q^{e',E)Fd{E)+ QBie',E)F^{E)]dE, (5) 

where we use the same cross-sections as in Equa- 
tion Q. 

The total observed flux is the sum of direct and 
back-scattered X-rays, i.e., /(e) — Io{^) + ^r(e) 
with /o(e), Ir(e), given by Equations ([2]), ([5]). Even 
using a single spacecraft we thus always observe a 
combination of reflected downward flux and direct 
flux which are functionally very different. This 
fact, and the properties of bremsstrahlung emis- 
sion, allow us to draw, from spectral data, conclu- 
sions about both the directivity and the spectrum 
of the emitting electron distribution. 

Figure ^ shows that albedo sets a natural 
constraint on possible directivity present in solar 
flares. Strong downward directivity (all electrons 
are confined within a pitch angle of 45°) typical to 
the collisional models (e.g. Brown, 1972; MacKin- 
non and Craig, 1992) leads to flatter than observed 
spectral indices below 20 keV and softer than ob- 
served spectra in the hundreds of keV range. It 
can be seen that the observed spectral index at 
200-300 keV cannot be less than 5 (Figure H]), 
much softer than typical spectral indices 2—3 (Ves- 
trand et al, 1987; Li et al, 1994). 



3. Tw^o directional electron distributions 

In reality we deal with discrete data sets rather 
than continuous functions. Therefore, the data 
vector of photons I(ei) = Iq + Ir can be written 

I = ( + G{^,)q^ + G{y)qr 

"(6) 

where e,; for i — 1...N, and Ej for j = 1...M, 
where Fd,u{Ej) are the electron data vectors and 
G, Q^'^ are matrix representations of the ker- 
nels of integral equations ()2I5|) . Green's matrixes 
G(/i) depend on heliocentric angle of the source 
fi = cos{9) and have been calculated in Kontar et 
al, 2006. The multiple sum in Equation © E*^E2 
can be expressed as a single sum S^*^ (Hubeny 
and Judge 1995). Equation ([6]) can be solved us- 
ing the Tikhonov (1963) regularization method of 
which our implementation was successfully tested 
by simulation (Kontar et al. 2004) and applied to 
RHESSI data (Kontar et al. 2005). 

To test our method we applied it to simulated 
data from known electron distributions, adding 
typical noise. Two different forms of anisotropy: 
'weak' anisotropy F„(£;) = _FdiE)/{l+y/{E ~ 10) /50) 
and 'strong' anisotropy Fu{E) = Fd{E)/{l -\- 
y^{E — 10)/50)2 in downward direction. From 
the simulated photon spectra we inferred regular- 
ized electron fluxes {Fu{E) ,Fd{E)) . The results 
demonstrate that the method used recovers reli- 
ably the directional electron spectra. However, the 
method is not sensitive to too 'weak' anisotropy 
that gives us anisotropy sensitivity. The difference 
between Fu{E) and Fd{E) for the case of 'weak' 
anisotropy is within the error bars. 

First we consider some limiting cases. 

If all electrons were purely downward di- 
rected (e.g. Fd ^ Fu — in Equation ^) 
then the observed spectrum would be mostly given 
by reflected and backward emitted photons (Fig- 
ure 2). Since the efficiency of Compton scatter- 
ing as well as the efficiency of backward emis- 
sion {QB{E,e) « QFiE,e) for E,e > 50 keV) 
decrease fast with energy, the photon spectrum 
should be rather steep. The photon spectral index 
should be not less than 5 (Figure 1) in the 200-300 
keV energy range even for a flat power-law electron 
spectrum while the typical spectral index observed 
is around three (Aschwanden, 2002). Moreover, 
the inverse approach concludes that for the Au- 
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Fig. 1. — Simulated photon flux spectra (upper 
panel) and spectral index (lower panel) for a flare 
located at heliocentric angle 6 = 40° with elec- 
tron distribution {Fd{E) ~ E-'^-^ and Fu{E) = 0, 
e.g., all electrons are directed downward and con- 
fined within pitch angle of ±45°). Upper panel: 
Observed (solid line), downward directed (dotted 
line), upward directed (dash line) and reflected 
(dash-dotted line) flux spectra. Lower panel: ob- 
served (solid), upward (dotted), and downward 
(dashed) spectral indexes. 
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Fig. 2, — Lower panel: The recovered mean elec- 
tron flux spectra (thick lines) for August 20, 2002 
flare (accumulation time interval 08:25:20-08:25:40 
UT) downward-directed Fd{E) (solid line) and 
observer-directed Fu{E) (dash line) with corre- 
sponding errors (thin lines). Upper panel: elec- 
tron anisotropy defined as Fd{E)/Fu{E) with con- 
fidence values within shaded area. 
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gust 20, 2002 flare (Kasparova et al, 2005), if the 
HXRs were strongly downward beamed (within a 
pitch angle of 45°) the electron spectrum required 
would have to grow above 400 keV to fit the data 
which casts doubt on the original assumption of 
downward beaming. 

If we ignore albedo Equation ^ leads to 
a substantial gap or low energy cutoff in the 
mean electron spectrum (e.g. Kasparova et al, 
2005), hence the local speed distribution function 
f{v) = F{E)dE/dv/v has a positive derivative. 
Such distributions are unstable and relax in the 
solar corona on a timcscalc much shorter than our 
observing time interval (e.g., Emslie and Smith 
1984, Melnik et al, 1999). 

We now consider the spectra at the peak of two 
solar flares: August 20, 2002 around 08 : 20 UT 
(heliocentric angle ~ 40°) and January 17, 2005 
around 09 : 40 UT (hehocentric angle - 33°). 
Both have sufficiently high count rates for good 
statistics up to above 500 keV but not so high as 
to cause severe pulse pile-up effects (Smith et al, 
2002). We used 7 out of 9 front RHESSI segments, 
excluding detectors 2 and 7 due to their low en- 
ergy resolution at the time of observation (Smith 
et al, 2002). Naturally we have chosen events not 
far from the disk centre since limb events are not 
suitable for our analysis because albedo is almost 
negligible for such flares. 

The general properties of the August 20, 2002 
flare ~ 08 : 20UT, which matched our count rate 
criteria, were extensively analyzed in Kasparova 
et al, 2005. To fit the hard X-ray spectrum with 
a double power-law electron spectrum with a low 
energy cut-off Ec and ignoring albedo, requires an 
unusually high value of Ec ^ 30 ± 2 keV. This 
produces a clear gap in the overall F{E) in the 
range 15-30 keV which is likely to be unphysical 
and suggests albedo is important. 

Figure [2] shows electron spectra solutions 
{Fd{E) , Fu{E)) with the confidence strips of our 
inferred downward and upward electron flux spec- 
tra. The results are close to consistent with 
isotropy up to 100 keV with some indication 
of upward anisotropy above 100 keV. Purely 
isotropic distribution would be given by unity, 
e.g. FuiE) = Fd{E). The January 17, 2004 flare 
has similar characteristics (Figure [3]): a rather flat 
spectrum (photon spectral index < 2.3) and good 
count rate up to a few hundred keV. However, 



this flare started in the soft X-ray tail of another 
event. To avoid difficulties with background sub- 
straction, we used data only from 15 keV. 

4. Discussions and conclusions 

The analysis of RHESSI spectra shows that we 
can infer simultaneously information on the direc- 
tional and the energy distributions of electrons. It 
is clear that for the January 17 flare (Figure [3|) 
that the Fd{E), Fu{E) solution are so overlapping 
and ratio Fd{E) / Fu{E) is so close to unity that the 
electron distribution is consistent with isotropy at 
all E. For the August 20 Hare in Figure O the 
distribution is consistent with isotropy at low E 
(below 100 keV) but may be slightly beamed to- 
ward the observer at high E. Since our recov- 
eries are means about just 2 directions, a wide 
range of actual /i distributions are consistent with 
them. For example, the August 20 event some- 
what resembles a 'pancake' distribution of parti- 
cles circling perpendicular to the magnetic field, so 
that perpendicular component of electron energy 
should be slightly larger than parallel one (Figure 
[2]). However, none of these results resembles the 
beam-like form expected for the basic CTT model. 
For example, a simplified mean particle treatment 
(Brown, 1972) has no electrons propagating up- 
ward, while more detailed dispersive treatments 
(MacKinnon and Craig 1991) suggests downward 
anisotropics of up to Fd{E)/Fu{E) ~ 10 - 100 de- 
pending on energy and on initial pitch angle distri- 
bution. Full Fokker-Planck models including mag- 
netic mirroring with a converging magnetic field 
(Leach and Petrosian 1981) give larger upward 
flux but still predominantly downward and out- 
side our confidence interval. For the large fluxes 
involved, these events may involve substantial re- 
turn current _E-flelds which can produce some up- 
ward 'beaming' but the effect is larger at lower be- 
low 70-75 keV, not higher, energies (Zharkova and 
Gordovskyy 2005). More advanced self-consistent 
models with collective effects included are needed 
for detailed comparisons with the observations. 

In conclusion, when allowance is made for the 
albedo contribution to the observed HXR spectra, 
the absence of a strong albedo feature precludes 
the basic model with strong downward beaming 
including CTT models, at least for the two intense 
hard solar disk events we have analyzed. This 
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casts serious doubt on the CTT model, at least 
in terms of the usually adopted geometry and on 
its physical basis. Alternatives to the CTT Stan- 
dard Model may need to be considered, in which 
acceleration occurs in such a way that collimated 
injection into a cold target is not involved. For 
example, HXR source electrons could be locally 
and continuously reaccelerated with near isotropy 
along the entire length of a magnetic loop. This 
model may also solve some of the other difficul- 
ties associated with the CTT model such as the 
problematic high beam density involved. 
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Fig. 3. — The same as Figure ([2]) but for Jan- 
uary 17, 2005 flare (accumulation time interval 
09:43:24-09:44:20 UT). 
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